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ABSTRACT 

In this paper, we have discussed logistic mathematical model, to analyze the population problems of Bangladesh. 
We have also discussed analytic solution, geometrical analysis, stability and local linearization of the logistic growth 
model. We analyzed the logistic model, both analytically and numerically. Several census data are analyzed for some 
countries, including Bangladesh. We have analyzed the population growth of Bangladesh, for a certain period of time 
(1971-2011). Finally, we implement the logistic model, that gives the future population of Bangladesh during (2012-2050). 

KEYWORDS: Analytic Solution, Carrying Capacity, Equilibrium Point, Geometrical Analysis, Logistic Population 
Model, Local Linearization, Stability 

1.0 INTRODUCTION 

The logistic model was named in 1844-1845, by Pierre Francois Verhulst, who studied it in relation to population 
growth [1 l Verhulst derived his logistic equation, to describe the self-limiting growth of a biological population. 
The equation was rediscovered in 1911 by A.G. Mckendrick, for the growth of bacteria in both and experimentally tested 
using a technique, for nonlinear parameter estimates [2 l The model is also sometimes called Pearl, R., & Reed, L. J. 
(1920), on the rate of growth of the population of the United States, since 1790 [31 . The carrying capacity k (l ) leads to a 
logistic delay equation, which has a very rich behavior 141 . The main point about the logistic model is a particularly 
convenient form to take, when seeking qualitative dynamic behavior in populations in which N = 0 is an unstable steady 
state and N (t) tends to a finite positive stable steady state 151 . The logistic form will occur in a variety of different 
contexts, throughout the paper primarily, because of its algebraic simplicity and because, it provides a preliminary 
qualitative idea of what can occur, with more realistic forms. It is instructive to try to understand, why the logistic form 
was accepted since, it highlights an important point in modeling in the biomedical sciences. The logistic growth form in 

this paper contains three parameters, N 0 ,kand r with which, to assign to compare with actual data. 

These were used by Pearl (1925), to fit the census population data, for various countries, including the United States and 
France for various periods. The main point is not that, the predictions are so inaccurate, but rather that curve fitting, only of 
the data, and particularly the part which does not cover the major part of the growth curve makes comparison, with data 
and future predictions extremely unreliable. The population problem is one of the main problems in Bangladesh at the 
current time. Bangladesh is an overpopulated country and the growth in resources has not been keeping pace with the 
growth in population. So, increasing the trend in population is a great threat to the nation. Kabir and Chowdhury 
investigated the relationship between population growth and food production in Bangladesh [6]. Recognizing the difficulty 
of feeding the growing population, even with considerable increase in food production, they suggested giving priority to 
population policy, for reduction in population. In this situation, prediction of population is very essential for planning. 
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2.0 Formulation of the Logistic Model 


Considering the fact that the growth rate actually depends on the population, we replace the constant 
- in dN _ by a function h(N) and thus, we obtained the modified equation dN 


dt 


dt 


- h(N)N 


( 2 . 1 ) 


We now choose h(N) so that h(N) = r > 0 when N is small, /?( N ) decreases as N grows larger, and h(N) < 0 when 

N is sufficiently large. Having this property, the simplest function is h(N) = r — aN , where tl is a positive constant. 

This (2.1) becomes df L = {r _ aN)N (2 2 ) 

dt v ' ’ 

which is known as Verhulst equation or the Logistic equation. Its solution is called the logistic function. The graph of a 

logistic function is called a logistic curve. 

[P.F. Verhulst (1804-1849) was a Belgian mathematician who introduced an equation (2.2) as a model for human 
population growth in 1838. He referred to it as a logistic growth; hence equation (2.2) is often called the logistic equation. 
Due to inadequate census data he was unable to test the accuracy of his model.] 

Suppose an environment is capable of sustaining no more than a fixed number k 0 f individuals in its population, 
the quantity is called the carrying capacity of the environment. 

Hence the function h in (2.1) we have h{k)= 0, and we simply let h(o)= r ■ The simplest assumption that we 
can make is that h(N) is linear, that is, h(N) =c { N +c 2 . Then h(0) = r givesCj.O+c 2 = r =>c 2 = r. 


And h(k)= Ogives c k +c 0 => _£jL = _L 
12 k k 

-r N 

.. h(N) = -A + r = r(l-) 

k K 


(2.3) 


Where 


The equivalent form of the logistic equation is-= r( 1- )N 

dt K 

_ i_ The constant r is called the intrinsic growth rate, that is the growth rate in the absence of any limiting point. 

a 


3.0 Carrying Capacity 

The carrying capacity of a biological species in an environment is the maximum population size of the species 
that the environment can sustain indefinitely, given food, habitat, water, and other necessities available in the environment. 

3.1 Analytic Solution of Logistic Model 

Verhulst suggested the following model. 


dN „ N ^r 

— =r(\-—) N 

dt K 


(3.1) 


Where N = N (t) = population of a single species at time t, r and k are positive constants. 
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Here r n—)is the per capita birth rate. The constant k is the carrying capacity of the environment. From (2.3) 
k 


we have 


dN 

N 0 - ^) 
k 


= rdt 


kdN 


N(k - N ) 
1 1 


= rdt 


dN = rdt 


N k-N , 

f— dN + f dN = f rdt- [By integrating] 

1 N 1 k-N 1 

> In N — ln(ft - N) + In A = rf;ln A is integrating constant. 

, AN 

In-= rt 

k-N 

AN = e « 


K-N 

=>AN=(k—N)e" 

Initially at t = 0 if n=N(o)=N 0 then from (3.2) we get, 

AN 0 -(k-N 0 )e° 

. k-N 0 

A =- - 

N 0 

By putting this value in (3.2) we get, 
k ~ Nr '-N = {K-N)e rl 


(3.2) 


N n 


kN 


-N + Ne r ‘=ke rl 


> kN-NN 0 +NN 0 e r, = kN 0 e rt 


>N[k+N^e n -l)] = kN 0 e n 


N = 


kN 0 e r 


k +N {) (e n -\) 


N(t) = 


kN 0 e n 

k+N 0 (e n -l) 


(3.3) 
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Corollary-1: When t —^ 00 then N (?) -A k 
Proof (Corollary-1): From the equation (3.3) we have 

kN n e n 


N(t) = 


k+N 0 (e n - 1) 


kN n 


=> (V(f) = - 

ke~ r, +N 0 ( l-e~ n ) 

When t —^ 00 then N(t)^> 


kN n 


k.0+N 0 {l-0) 


■ N(f )->^ 


iV(0 —> k 


That is, we have shown that when time is unlimited, then the population approaches to the maximum, 
the equilibrium number of organisms of a particular species that can be supported indefinitely in a given environment. 


Corollary-2: Assume that are equally spaced time values. Suppose the corresponding population size are 

N q, and N , respectively, than the carrying capacity kof the population that obeys the logistic law of population 


can be expressed as N 0 N[ + N[N 2 2N 0 N t N 2 

N?-N 0 N 2 

Proof (Corollary-2): From the equation (3.3) we have AN=(k—N)e n (3.4) 

Initially ^^o’^t’^2 anc * N=N 0 ,N 1 ,N 2 respectively. So we have from (3.4) 

AN 0 =(k-N 0 )e rt ° (3.5) 

AN l =(k-N l )e rt ' (3.6) 

AK=(k-N 2 )e h (3.7) 

(3.5) - 7 - (3.6) gives ^0 = k ~ N o e ^-h) (3.8) 

N { k-N l 

(3.6) A (3.7) gives N ' = k ~ N t e r ^ (3.9) 

N 2 k—N 2 


Since ^o , ^l’^2 are equally spaced, so we have t ff ^ 
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Thus, (3.9) becomes = ^L^JL e 'i>o h) (3.10) 

(Vj k-N 1 

(3.8) H- (3.9) gives, N o N 2= k-N 0 k-N 2 
(V, N l k-N l k—N l 

_ N 0 N 2 _ (k-N 0 )(k-N 2 ) 

Nf ( k-Nf 


=>N 0 N 2 (k -Nf=Ni 2 (k -N ( ) (k -N 2 ) 
^>(k 2 -7kN x +Ni 2 )N q N 2 =Ni 2 (k 2 -kN-kN Q +NJV 2 ) 


=>k 2 N 0 N-2kN 0 N 1 N 2 +Ni 2 N 0 N 2 =k 2 Ni 2 -kNx 2 N 2 -kN 0 Ni 2 +N 0 M 2 N 2 
=>.k\N 0 N 2 -Ni 2 )= k(2N Q (V, N-Ni 2 N-NJf i 2 ) 

=>k(N 0 N 2 -Ni 2 )=2N 0 N 1 (V 2 -(Vi 2 A^ 2 -A^oA^i 2 [v k* 0] 

K = N 0 N 2 +N 2 N 2 -2N 0 N l N 2 

N i~ N o N 2 (Proved). (3.11) 

(j__ 0 

Corollary-3: Suppose q < A, then j N £ I 

0 2 r = -l n 7 — -f 

f n _i 

U. 

Proof (Corollary-3): From the equation (3.3) we have N( ., kN 0 e" 

' k+N 0 (e r '-l) 

At time ^, let )V(f) = 2(V 0 


A(f) = 


kN „e" 


k+N 0 (e*- 1) 


- = 21V„ 


=>ke r '=2k+2N 0 e r, -2N 0 

MK-2N 0 )e r '=2{k-N 0 ) 

n _ 2(k-N 0 ) 

k-2N 0 

, 2(*-JV 0 ) 

=> rt = In—-— 

k-2N 0 

1 , 2(k-N 0 ) 

=>t = - In—-°' 

r £-2)V 0 
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r = 

t k-2N 0 


=> r = 



t 


2(k-N 0 ) 



=> r = 



t 


{k ~N o) 



Suppose & = — =>#=! 

N o * 

Let us consider 2N 0 =N ] 


1 _J. 

.,.' 4^1 

f M_J. 

u * 


4.0 Obtaining Equilibrium Points 

We obtain the system’s equilibrium points iV * by finding all values of N that satisfy _ q . 

dt 


0 =>^r*ft—^3=0 (4 - 1} 

dt K 


.-. rN * =0 or| 1-— | = 0 
AT 


We get N*= 0 (4.2) 


N*= K (4.3) 


Thus, the logistic equation has exactly two equilibrium points. 

Table 4.1: Behavior of Logistic Growth, Described by Equation (3.1) for Different Cases of N 


N 

dN 

dt 

N> K 

0 

dt 

0<N <K 

">0 

dt 

N = K 

™=0 

dt 

N=0 

dt 
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5.0 Geometrical Analysis 

Table 4.1 illustrates the sign of for different values N . The trivial equilibrium point N * is unstable, and 

dt 

the second equilibrium point N = k represents the stable equilibrium, where N asymptotically approaches the carrying 

capacity K . In terms of the limit, we can say ii m jV(T) = k N( 0) > 0 A point of inflection occurs at _ k for all solutions 

T^°° ~2 

that cross it, and we can see graphically that.owth of N is rapid until it passes the inflection point _ k . 

~ 2 

From there, subsequent growth slows as N asymptotically approaches k, as shown in table 4.1, HN>k, then dN 

dt 

and iV decreases exponentially towards k. This case should only occur when the initial condition A^o^lV^A:. 
In the following section, we will confirm the stability of equilibrium by linearization about each equilibrium solution. 


dN 

dt 



Figure 5.1: Phase Line Portrait of Logistic Growth, as Described by Equation (2.3) 


N 



Figure 5.2: Dynamics of Logistic Growth, as Described by Equation (3.3) 
6.0 Local Linearization 

Let us consider a population in general to be governed by 


dN 

dt 


f(N ) 


( 6 . 1 ) 


Where f(N)is a continuously differentiable (typically nonlinear) function ol'/V then the equilibrium 
solutions ./V * are solutions of f{N) = 0 and linearly stable to small perturbations if /'(a^‘)<0 . unstable if f'[N*)> 0. 
This is clear from line arising about N*, by writing n(t) ~ Nit) —N* , n(t )\«I • Then (6.1), 

becomes cl *L= f( N *+n)» f(N*) + n f\N*) +... which to first order in n(f)gives ) =><0 °= enp[fXN*)t] (6.2) 


dt 


dt 


So kl grows or decays accordingly as /'(Al')>0or / / (w)<0- Timescale of the response of the population to a disturbance is 
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of the order of 1 It is the time to change the initial disturbance by a factor C. Depending on the system ./ {N ) 

j/'H’ 

models, may get several equilibrium or steady state populations N * . Graph plotting f (/V ) against N immediately gives 
the equilibria as the points, where it crosses the N -axis. The gradient f\N') at steady state, then determines its linear 
stability. Such steady states may, are positive, so these equilibria are unstable are stable to small perturbations: the arrows 
symbolically indicate stability or instability, so that, N is in the range /V,< N <N, then N —>/V, rather than returning to 

N j. A similar perturbation from to a value in the range 0 < N <N , would result in N(t)=>N r 

Qualitatively, there is a threshold perturbation below, which the steady states are always stable, and this threshold depends 
on the full nonlinear form of f(N). For example, the necessary threshold perturbation is N 2 -N r 


N 

Figure 6.1: Population Dynamics Model 4K _ j iN j with Several Steady States. The Gradient f'(N) at the Steady 

dt 

State. That is, where f(N) = 0, Determines the Linear Stability 
7.0 Realistic Analysis of Logistic Model Used to Fit the Census Data 

In this paper, we have discussed about the logistic model by Verhulst (1838) with its realistic situation in some 
countries. We have considered here the countries Belgium, USA, France and Bangladeshi populations for a certain period 
of times and analyzed their numerical solutions with the model. We have shown here that this model fits for either lower or 
upper parts of the graph-model for some countries whereas all of the lower or upper part fits for a certain country like 
Bangladesh and USA. 

7.0 Logistic Model Used to Fit the Census Data for Belgium 

Using the Corollary (2) and corollary (3) estimation of the population of Belgium in the years 1815, 1830 and 
1845 (respectively 3.627, 4.247 and 4.801 million), we obtained K =6.584 million and r = 2.62% per year. 
The population of Belgium and the logistic curve showing the following Figure (7.1). 



Figure 7.1: The Population of Belgium (in Million) and Logistic Curve 
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7.2 Logistic Model Used to Fit the Census Data for USA 

Pearl (1925) used logistic model to fit the census data of USA population for various periods. 
Figure 7.2 shows the results for USA. If we look at the USA data there is a good fit for the population roughly from 1790 
until about 1910; here the lower part of the curve is fitted and later we shown a table (Table 7.1) of USA census data and 
analytic value in the period (1970-2010). 



Figure 7.2: Logistic Population Growth Eq. Used to Fit the Census Data for the Population of USA (In Million) 


Table 7.1: USA Census Data and Analytic Value Since (1790-2010) 


Year 

USA Population 

Analytic 

1790 

3,929,214 

3,929,214 

1800 

5,236,631 

5,335,741 

1810 

7,239,881 

7,227,100 

1820 

9,638,453 

9,755,109 

1830 

12 ,866,020 

13,106,950 

1840 

17,069,453 

17,503,927 

1850 

23,191,876 

23,191,876 

1860 

31,443,321 

30,418,277 

1870 

38,558,371 

39,391,703 

1880 

49,371,340 

50,223,418 

1890 

62,979,766 

62,860,165 

1900 

76,212,168 

77,029,658 

1910 

92,228,531 

92,228,530 

1920 

106,021,568 

106,124,688 

1930 

123,202,660 

121,713,981 

1940 

32,165,129 

39,081,412 

1950 

151,325,798 

158,280,419 

1960 

179,323,175 

179,323,175 

1970 

203,211,926 

202,171,349 

1980 

226,545,805 

226,728,483 

1990 

248,709,873 

252,835,188 

2000 

281,421,906 

280,268,397 

2010 

308,745,531 

308,745,531 


7.3 Logistic Model Used to Fit the Census Data for France 

Verhulst in 1938 did a similar study of France population, finally Pearl in 1925 used logistic model, to fit the 
census data of France population for various periods. Figure 7.3, shows the results for France. If we look at the France 
data, there is only fitted to the upper part of the curve, but even there the subsequent population growth prediction is 
wrong. 
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year 

Figure 7.3: Logistic Population Growth Eq. Used to Fit the Census Data for the Population of France (In Millions) 
7.4 Logistic Model Used to Fit the Census Data for Bangladesh 


To fit the census data of the Bangladesh population, we need the value of carrying capacity K and growth rate r . 

We have used corollary (2) and corollary (3) to find carrying capacity K and growth rate r of Bangladesh population. 
Using Verhulst procedure that means in the help of corollary (2) and corollary (3) we show table (7.2) Bangladesh actual 
population, analytic and numeric data using Runge-Kutta scheme. 

Table 7.2: Census, Numeric and Analytic Data of Bangladesh Population 


Year 

Bangladesh Population 

Numeric 

Analytic 

1971 

7.16 



1976 


8.083 

8.083 

1981 

8.98 

9.069 

9.069 

1986 


10.095 

10.095 

1991 

11.16 

11.145 

12.203 

2001 

13.06 

13.250 

13.250 

2006 


14.271 

14.272 

2011 

15.26 

15.252 

15.251 


(Population in million) 



Figure 7.5: Logistic Population Growth Equation Used to Fit the Census Data for the Population of Bangladesh 

(In Millions) 
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Figure 7.6: Logistic Population Growth Equation Used to Fit the Census Data for the Population of Bangladesh 

(In Million) by Pie Diagram 

8.0 Future Population of Bangladesh Using Logistic Model 

Using logistic model (2.3) we show table (8.1) Bangladesh population data from 2012 to 2050 is increasing 
without bound. 


Table 8.1: For the population of Bangladesh (2012-2050) 


Year 

Population of Non Linear Model 

2012 

15.0425 

2013 

15.2501 

2014 

15.4611 

2015 

15.6761 

2016 

15.8960 

2017 

16.122 

2018 

16.3556 

2019 

16.5986 

2020 

16.8531 

2021 

17.1215 

2022 

17.4065 

2023 

17.7111 

2024 

18.0386 

2025 

18.3921 

2026 

18.7745 

2027 

19.1885 

2028 

19.6366 

2029 

20.1211 

2030 

20.6441 

2031 

21.2075 

2032 

21.8130 

2033 

22.4621 

2034 

23.1561 

2035 

23.8961 

2036 

24.6830 
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Table 8.1 

Year 

Population of Non Linear Model 

2037 

25.5175 

2038 

26.4001 

2039 

27.3311 

2040 

28.3106 

2041 

29.3385 

2042 

30.4145 

2043 

31.5381 

2044 

32.7086 

2045 

33.9251 

2046 

35.1865 

2047 

36.4915 

2048 

37.8386 

2049 

39.2261 

2050 

40.6521 


(Population in million) 


pcpJaticn 


Vears 

2012 2016 2020 2C24 2028 2032 2036 2040 2014 2)18 2050 

Figure 8.2: Logistic Population Growth Equation Used to Fit the Census Data for the Population of Bangladesh 

(In Millions) 



Figure 8.3: Logistic Population Growth Equation Used to Fit the Census Data for the Population Bangladesh 

(In Million) by Pie Diagram 

9.0 CONCLUSIONS 

We have discussed in this paper of some basic concepts, dynamical behavior and numerical solution of the 
logistic mathematical model based on population dynamics easily. We analyzed the data of Belgium, USA and France. 
Then we used the same method to project population for the aforementioned period of Bangladesh. We observed that our 
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data for the same period of time fitted with a good agreement with the census data. In this study a mathematical analysis of 
the future population of Bangladesh is carried out based on an ordinary differential equation model which is called logistic 
model. MATLAB program based on an algorithm of Runge-Kutta scheme is developed for direct calculation of future 
populations. Then we establish a non-linear model that gave future population in Bangladesh at the time from 1971 to 2011 
and 2012 to 2050. The realistic mathematical logistic model is useful in many applications from experimental instruments 
to rigorous mathematical analysis techniques. Finally, we have noticed that the population of Bangladesh is increasing 
without bound. 
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